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Although combat is a complex random process, analysts 
frequently model it with deterministic Lanchester-type 
equations for reasons of mathematical tractability and com- 
putational convenience. Such deterministic models are being 
used extensively today in defense planning studies. However, 
one can include the various random aspects of combat (e.g., 
uncertainty in enemy force level, random occurrence of 
casualties, etc.) in models and develop so-called stochastic 
models of the combat process. Thus, a basic way of classi- 
fying a combat model is whether the model is a deterministic 
model or a stochastic one. 

Moreover, there are fundamental differences between 
deterministic and stochastic models. Deterministic dynamic 
combat models predict the future with certainty; for given 
initial conditions there is no question about what the state 
Seetne conflict will be at aay, future time. On the other 
hand, stochastic dynamic models only tell one the "chances" 
of what will happen. With a stochastic model one does not 
know the future with certainty. That is, it cannot be 
guaranteed that a specific state will be reached in the 
future — there is only a "probability" it will be reached. 

Each type of model has inherent advantages and disadvantages 


which are partially summarized in Table I. 
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Deterministic Models 


Advantages Disadvantages 
* 1. Requires one solution of the IE ae (ire [cl Veh aera one mrerclehaulre}el 
model for a useful outcome from reality 


2. Less expensive to run 


3. Interactions are more 
easily analyzed 


Stochastic Models 


Advantages Disadvantages 
1. Closer to reality 1. Requires many replica- 
tions to provide 
statistically 


significant results 
2. More expensive to run 
3. Interactions are much 


mMOLre CYrereult tO 
analyze 


TABLE I. Advantages and Disadvantages 
of Deterministic and Stochastic 
Models 
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The solutions tostochastic models quickly become quite com- 
plex as the number of combatants increases, and are not very 
useful; so one must use other methods to obtain insights into 
the dynamics of combat. The two most commonly used methods, 
Monte Carlo simulation and finite difference approximation, 
both possess certain disadvantages. Monte Carlo combat simu- 
lations require that the battle be replicated many times in 
order to obtain a good statistical estimate of the probable 
course of battle, while finite different approximation methods 
usually require digital computer implementation (including 
development of a computer program). 

On the other hand, solutions to deterministic models can 
sometimes be useful for developing insights into the dynamics 
of combat. Even when they are not, the model must be solved 
only once by finite difference approximation as opposed to 
the many replications required for Monte Carlo simulation. 
Also, finite difference approximations of deterministic models 
are generally less costly in terms of time, and manhours 
required to implement the computer program, than the finite 
difference approximations of stochastic models. The economic 
costs of using the models then, is the basic reason most models 
of combat used in the past have been deterministic. But can 
the complex random process of combat be adequately represented 
by a deterministic model? Does one destroy any of the essential 
features of the random combat process by considering a deter- 


ministic model as representing the "average" course of combat? 
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The objective of this thesis is to consider this question by 
investigating the differences in results obtained from a 
deterministic and a stochastic model of combat. 

The importance of this objective stems from the fact that 
both types of models are used extensively today in the decision 
making process in the U.S. Army. For example, the Dynamic 
Tactical Simulation, DYNTACS-X, and the Battalion Differential 
Model, BLDM, widely known as the Bonder/IUA model, are two 
mode 1s that simulate battalion level mid-intensity armored 
combat. DYNTACS is an event Sequenced stochastic model that 
memoerng Or has been used by the Rock Island Arsenal to evalu- 
ate proposed mobility improvements to the M60 tank, the family 
of scatterable mines concept, cannon launched guided projec- 
tiles, remotely piloted vehicles, and the XMl tank. Bonder/ 
IUA is a deterministic Lanchester-type model that is being or 
has been used by the Rock Island Arsenal in sensitivity analyses 
and evaluation of the MBT70 study, several anti-armor automatic 
cannon concepts, and the Low Dispersion Automatic Cannon study. 
A study done by the Rock Island Arsenal [Ref. 5] shows that 
there are 4s Eons in the results of the two models with 
"equivalent" inputs. With the present state of the art, it 
is not possible to predict with any accuracy when significant 
differences between the results of the models will occur. It 
is felt that this thesis may provide trends which will allow 
a user to gain insights into these differences and more 
accurately determine when a deterministic model represents 


Someat “Satisfactorily". 
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Models have had wide application in the military decision 
Making process, and the realm of possible future applications 
is even larger. For example, the Engineer Strategic Studies 
Group [Ref. 10] has, as itsS name implies, evaluated several 
aspects of military engineering, but with little emphasis on 
the use of combat models. Much of this has been caused by a 
dearth of models that could be applied to engineer studies. 
However, this is being altered; for example, the U.S. Army 
Concepts Analysis Agency [Ref. 30] has accepted as a future 
tool a model developed by ESSG, the Fortification and Obstacle 
Mecect Simulation, 1975. This is a high resolution model 
used in analyzing obstacle effectiveness. The model is part 
of a research effort by ESSG to design a method for assessing 
the contribution that obstacles make to combat operations. 
This model was developed using ideas from other models, but 
1t was based on the identification of specific objectives. 
From these objectives measures of effectiveness were developed, 
and from these measures the model was developed. It is possi- 
ble, with the identification of proper objectives, to develop 
models for other phases of military engineering and, more 
generally, for the evaluation of other phases of combat. It 
is therefore felt that combat models do, and will continue to 
have an ever increasing application to the evaluation of 
Military systems. 

As was noted, however, it is not always clear which type 
of model should be applied — a stochastic or a deterministic 


model. From an economic standpoint, the advantages of the 
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deterministic model weight heavily in its favor. If the 
deterministic results do not significantly differ from those 
of the stochastic model, it seems the former should be used. 


The question is, when are the models "Similar"? 


Li 





II. LITERATURE SURVEY 


This chapter briefly reviews Lanchester's well-known 
equations of combat and some simple extensions, with the 
discussion focused on stochastic formulations of the combat 
attrition process. The review consists of three parts; 

(A) Deterministic Models, (B) Stochastic Models, and 


(C) Comparisons of the two types. 


a DETERMINISTIC MODELS 
In 1914 Lanchester [Ref. 19] hypothesized that under 
conditions of "modern warfare" combat between two homogeneous 


forces could be described by the equations 


cos 
ale es = 
DY ae 
Aig bx (2) 


with initial conditions 
x(t = 0) = x ¥35) 
y(t = 0) =y, (4) 
where a and b are commonly referred to as the Lanchester 
tration rate coefficients and x(t) and y(t) are force 


levels. One set of assumptions that has been hypothesized 
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(proposed by H. Weiss [Ref. 35] and cited by Dolansky 
a 
{[Ref. 9]) to yield equations (1) and (2) is 
1. Two forces are engaged in a battle. The units 
are homogeneous, but the rate of attrition may 


be different for each force. 


2. Each firer on either side is within the range 
of all targets on the other side. 


See ttLition rate coefficients are constant. 
4(a). Firers have perfect knowledge of the target 
locations and disposition so that they fire only 


at live targets, and fire can be shifted 
instantaneously. 


5(a). Fire is uniformly distributed over remaining 
targets. 


The state solution (time "independent" solution relating 


the force levels ) is given by 


b(x,* - x(t)7) = aly.* - y(t)*) (5) 


hence the name "Square-law" attrition process. 
The time solution for equations (1) and (2) is given 


by (see Morse and Kimball [Ref. 22]}) 


x(t) = x, cosh (yj ab t) - 1 a/b y, sinh (\fab t) (6) 
We) = yeecosh (Gyabrt) - 1, b/a x, sinh (\jab t) (7) 


The square-law was developed as a by-product of the 
Original work by Lanchester to quantitatively justify the 


principle of concentration under certain combat conditions. 
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Lanchester was contrasting these conditions with the condi- 
tions of "ancient warfare" between two homogeneous forces. 
He postulated that the latter could be described by the 


equations 


<> a 
se ES ze 
chy aan 
ae bxy (9) 


Pee initial conditions 


x(t = 0) = x (10) 


y(t = 0) (27) 


i! 
he 


One set of assumptions that has been hypothesized (also 
proposed by H. Weiss [Ref. 6) and cited by Dolansky [Ref. 9]) 
to yield equations (10) and (11) are very similar to the 
assumptions listed previously for the square-law. The first 
three assumptions are identical to the first three square- 
law assumptions. 

4(b). Each firer knows the area that contains 

targets, but does not know exact target 


locations or the consequences of his fire. 


io). Fire £Erom Surviving units is distributed 
uniformly over the target area. 


The state solution is given by 


b(x, - x(t)) = aly, - y(t)) (12) 
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hence the name "linear-law" attrition process. 
The time solution for equations (8) and (9) is given 


by (see G. Weiss [Ref. 35]) 


(y ~b/ax,) exp( (b/ax-y_) at) 


x(t) = Wy =b/ax, exp((b/ax,-y,)at) wal 
y(t) = Y¥,7b/ax, + b/ax(t) (14) 


Others have suggested general forms of homogeneous force 


models given by 


St = ~- F(t,x,y) (15) 
gy = -—- G(t,x,y) (16) 


For example, Willard [Ref. 37] considered the following 


Lanchester-type equations: 


d c | 
a =- ax” £(t,x,y) (17) 
Gly e C 


Many "general" forms contain variable coefficients (e.g., 
a= a(t) and b = b(t)). Taylor [Ref. 30] has shown that 
except for some very special cases, the solutions to varia- 
ble coefficient Lanchester-type equations are very complex 


and in many cases of little practical use. 
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8. OfOCHASTIC MODELS 

Koopman [{Ref. 22] suggested a reformulation of the 
attrition process in stochastic form. The resulting sto- 
chastic attrition process has been appropriately called by 
him the Lanchester stochastic process. Others have subse- 
quently employed a stochastic analysis of combat. With the 
inclusion of random variations in the attrition process, 
analysts have attempted to better represent the complex 
random process of combat with the goal of gaining insights 
into combat not available from deterministic models. 

Snow [Ref. 27] shded that the following assumptions 


yield the square-law attrition process: 


anAt 


(Al) P(ttAt,m,n:t,mtl,n) 
oz P(tTAE,M,n:t,m,ntl) = bmAt 


(A3) Probability of more than one casualty 
GeeunrtIng 2m a time Of length At is of 
the order of magnitude o(h) where 

o(h) _ g 

1 : 





lim 
h>0 


From these, he developed the Chapman-Kolmogorov forward 


equations: 
cet emen) = anP(t,mtl,n) + bmP(t,m,ntl) - (an+tbm) P(t,m,n) 
(19) 


If one assumes, instead 


(Al) P(ttAt,m,n:t,mtl,n) = F(t,mtl,n)At 


eee (et At,m,n:t,m,nt+)) G Cee mya) At 
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(A3) Probability of more than one casualty 
Sceurring im a time Of length At is of 
Magnitude o(h). 


then the resulting Chapman-Kolmogorov forward equations are: 


adP(t,m,n) 


at = F(t,mtl,n)P(t,mtl,n) + G(t,m,n+l)P(t,m,ntl) 


Peet, i, n) + Gltem,n)) P(t, m,n) (20) 


In this thesis, this process will be called the general 
Memenester—type stochastic™attrition process. 

Work has been done by several authors to obtain general 
Pemuctions for specific stochastic attrition processes. 


Brown [Ref. 7] attempted to find the time state probabili- 


ties (P(t,m,n)) for a class of homogeneous force stationary 
Markov attritions processes (that is, when F(t,m,n) = A(m,n) 
and G(t,m,n) = B(m,n)). He concluded that the expression 


derived is too complicated for practical use. Isbell and 
Marlow [Ref. 14] developed a general solution to (20) for 

a square-law stochastic attrition process. If 

F(t,;m,n) = an + 8m and G(t,m,n) = bm + an (Isbell and Marlow 
referred to a and 8 as operational loss rates), with 


ata= b+ 8, then 


motno-m-n 
eet, m,n) = F(m,nim ,n imo tn \exp ( (b+Bt) -1) 


m+n (21) 


» exp (-(b+B) (m +n.) t) 
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where F(m,n:im ,n) is developed from a recursive relationship 


moran, 
F(m,n:m,n) = (reym Flatayn, FMM EM“ h No) 
ntbm, 
Se ee ° = 
F(m,nimo,n, 1) (22) 


(b+8)m +(ata)n. 


where F(mo eng ?mMo ny) = l. 
Clark [Ref. 8] developed a general solution for the time 


state probabilities for the linear-law stochastic attrition 


process (i.e., F(t,m,n) =amn and G(t,m,n) = bmn); 
Wl Se Si) 
Mo Po | (-1)8°F9"M ao pl (myr (ne)! 
eet 2) = re ; ( gy el il al 
Jam k=n (atb) mio nt (k-n)! (j-m)! (no-k)! 
| | 


\ 
a J | 
58 genet "8? ngs sas 
an (j-k- ) =i Tj+k+L) SES? Ve clan ey) IS 


(23) 


(However, he did not compare his results with those for the 
probability of winning (i.e., "true" absorption probabili- 
ties). This expression is too complicated to provide any 
insights by itself. This author has no knowledge of any 
other general solutions developed for the time state proba- 
bilities. It is, however, a straight forward task to compute 
P(t,;mo-n) and P(t,m,no) » which are used in the included 


Semputer program aS a partial check for accuracy. 
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Work has also been done to develop exact expressions for 
the probability of winning, for a stationary Markov attrition 
process. The exact expression for the probability of winning 
for the linear-law stochastic attrition process is easily 


obtained from classical random walk results [Ref. 28]: 


O O Tete lek) 1: O Oo 
. a oC 2 _D_ 
P= °F P(o,n:m ,n) = oN mo -1 ) an arm 


(24) 
Smith [Ref. 26] extended work done by Brown [Ref. 7] and 
developed the explicit expression for the square-law stochas- 


tic attrition process: 


m 
Oo 
Pg z P(m,O:m_ ,n,) (25) 
m=1 
where 
m-} m_+n_-m-1l 
P(m,Q:m Fag) ) = iia is pel) eS yee) 
S20 an +bm. j=m P(mo-j+1) P(n +b/a giaedl)) (| =iweale) 


(26) 
This is an impressive looking expression, but from a practical 
viewpoint, it is not very useful. A computer is needed to 
solve this equation for any appreciable force sizes, but 
One immediately runs into overflow and underflow problems. 
To solve these complications, recursive relationships must 
be used; but if one looks carefully at the relationships 
used, he is lead back to the initial recursive relationships 


used to obtain the solution: 


ZS 





A(m/n) 


EaCOe ni sn 0) 


a P(0,n:mo-1,n,) 


A(m>,n,) +B(mo,n,) 


B(m +n.) 


+ PiCO-n:m-,n.—)) (27) 
av a? 
A(m,N,) +B (m,n) o° oO 


It was found by this author that the most efficient method 
to solve the stochastic model for the probability of winning 
was by uSing the above recursive relationship, or one devel- 
oped by Springall [Ref. 28], not the general force level 
solutions (See Appendix A). 

Springall developed a recursive relationship.for the 


distribution of the time of battle termination: 


n 
Oo 
tr) (xr) 
E(t) = ie ee XS 41 
b Pp 
Pp 
m 
+ a B(m #1) x'*? (28) 
a 41 ‘bp m,n), +] 
"5p p 
where 
(0) 1 
x = AG Ce i) 
MoD, A(mo7N, +B(m,n,) 
(0) 
ry ot teen. ntl 
a n A(m_,n)+B(m_,n) (30) 
oO oy! row 
(0) 
A(m+1 No Xoa on 
4 7 4 
m,n, A(m,n,) +B (m a) 
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(0) (0) 
Atm+l,n) Xo, ytBimentl) x ay 


Xan Atm,n)+B(n) ia 


panne 
yx (2) _ oo’ oO (33) 
aD a a 
Moe, A(morng) tB(mo rn, 
».(i-1) (1) 
1X iS +B(mo-ntl) xX “piel 
x ft) S O O (34) 
moyen A(m_,n) +B (m,n) 
oe nel yx?) 
, os inte, i 
yx {2) LE O O (35) 
m,n, A(m,n,) +B (m,n_) 
. .., (1-1) (1) (1) 
pare = aa tA(m+1 n) Xo ptB(m ntl) x ay 
‘ A (m,n) ey pew, 
(36) 


From studying the equations presented in this section, one 
gains an appreciation for the increase in complexity of the 
model by the inclusion of uncertainty in the attrition 
process. It is apparent that, although exact solutions may 
be obtained for some stochastic models, numerical results 
are not readily generated by hand. A computer is therefore 
essential for the numerical solution of the simrle stochastic 
models. 

It is surprising that more use has not been made of the 


computer in investigating stochastic models, especially the 
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use of computer graphics. For example, the five following 
figures are plots of the time state probabilities for a 
stochastic model at 3, 25, 50, 75, and 100% of the time to 
completion of the equivalent deterministic model. It is 
believed that something of this nature has not been utilized 
in the open literature before. It gives some interesting 
insights into the progress of a stochastic model. The state 
probability seems to start as a spike at (morn) and pro- 
gress” over time somewhat as a drop of water would. As time 
increases, the drop spreads and moves toward the banks 
(boundaries of the state space). The angle at which the 
wave moves is determined by the relationship between the 
attrition rates of the process (in this case, the model is 
the square-law stochastic attrition process, so the rela- 
tionship is between a, b, m, andn). As the "wave" hits 

the wall, it "sticks" (is absorbed), thus giving, as t-*~, 
the distribution of win probabilities. (It is interesting 
me note that this relates the work of Clark [Ref. 8], 

Weiss [Ref. 35], Morse and Kimball [Ref. 22], and Snow 

[Ref. 27], who were concerned with the solution to the 
Chapman-Kolmogorov, to the work of Smith [Ref. 26], Brown 
[Ref. 7], etc., who were concerned with the win probabilities). 
Note that in the present case, Y has a very large probability 
of winning. If the forces were moved closer to parity, more 
of the wave would be absorbed by the Y force boundary, 


Giving X a larger probability of winning. 
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Plot of P(t,m,n) for Fixed t 
FICURE 5 





The above discussion has given a description of what 
occurs with the probability mass over time. Another des- 
cription of the process is, it acts as a combination of 
diffusion and convective transport of probability mass. 
This concept enhances ones understanding of the nature of 
so-called diffusion approximations of Markov processes 
(l.e., state space assumed continuous). Additionally, it 
provides insights into temporal changes in variance in force 
levels (see Chapter III, Section C). Thus, it is felt that 
Semputer graphics, if utilized more fully could provide 
eaSily obtainable insights into the dynamics of a combat 


moael!, that would be difficult to obtain otherwise. 


C. COMPARISONS 
Two types of comparisons have been made between stochas- 
tic and deterministic models; probability of winning, and 
the expected force levels against the time history of the 
force levels of the deterministic model. Two recent studies 
have been done comparing the probability of winning with 
the certainty of winning or losing in the deterministic 
model [{Ref. 20 and Ref. 28]. Both came up with basically 
the same conclusions. 
(1) As force levels increase, the probability 
of winning approaches the certain probability 
of winning or losing in the deterministic 


model. 


(2) Significant differences occur only when the 
Pmecc mate Neat Pari ty . 
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Lee and Wannasilpa [Ref. 20] came to their conclusions for 

the square-law, linear-law, and mixed-law. Springall arrived 

at the same conclusions for a slightly more complicated model. 

Snow [{Ref. 27] showed that the time history of the 

expected force levels differed from the time history of 

the deterministic force levels; that is m(t)-x(t) # 0. 

fer (Rell egiedetined this as bias. This seems to be an 

appropriate term and will be used here. Snow apparently 

first developed an explicit expression for this bias for 

the square-law; 
en =-anta pormeae (37) 

Clark [Ref. 8] empirically showed that bias exists for the 

linear-law. All researchers seem to have arrived at the 

Same general conclusions: 


(1) As initial force levels increase, the bias 
decreases. 


(2) The longer a battle lasts, the greater 
the bias. 


Powers and Taylor [Ref. 24] and Hanna [Ref. 12] have 
Mempered duftferences in optimal time sequential fire distri- 
bution policies for deterministic and stochastic models: 
Hanna's investigation indicates that there is a real differ- 
ence between the two (at least for small numbers of total 


combatants). 
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Thus, some work has been done to investigate the differ- 
ences between deterministic and stochastic models, especially 
in the areas of who is going to win and the average time 
histories of the force levels. However, the analysis done 
up to the present has been somewhat limited in scope and 
does not allow an analyst to say with any degree of cer- 
tainty, when the complex random process of combat can be 


adequately represented by a deterministic model. 
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Pi CcOUEmnRiooNe Or DETERMINISTIC AND STOCHASTIC 


MODE Eo eno SQUARE, LAW, ATTRITION PROCESS 


It would be most useful to be able to make comparisons 


between stochastic and deterministic models used today. 


This 1s not feasible within the scope of this paper for 


several reasons. 


le 


It is not clear, in many cases, when the two models 
are equivalent. For example, one could try and 
compare DYNTACS-X and Bonder/I.U.A., but before 

this comparison could be made, one would have to 
insure that the input parameters are equivalent. 

In the case of these two models, that is not a 
trivial task, as the structure of the two models 

is quite different, particularly in route selection. 
About the best that can be done is to construct the 
input parameters for each model and assume they 

are equivalent. 

Parametric analysis is not feasible in a complicated 
model. Complicated models, both stochastic and 
deterministic, are expensive to run. Because of 
this and the complex interactions among the variables 
of a model, parametric analysis is not only 
prohibitively expensive, but difficult to do. 

It is difficult to gain insights in a comparison 


of simple stochastic and deterministic models, and 
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essentially impossible to gain many useful insights 

in a comparison of complicated models. 
Lanchester-type models may be thought of as resulting from 
an aggregation of the myriad of combat environmental varia- 
bles into constant attrition rate coefficients and time 
varying force levels and, moreover, as complementary complex 
system models. Such idealizations will be considered in 
this thesis with the knowledge they are idealizations, but 
that they may provide insights that can be used as points 
of departure for the investigation of more complex models 
of combat. Specifically, two models will be compared; the 
deterministic Lanchester square-law and the square-law 
stochastic attrition process. It is believed that trends 
found from this investigation may be carefully used to 
determine trends for a comparison of more complex models. 

As was shown in the literature survey, the square-law 


1s a Special case of the general differential equations 


a Ge yt) (LS) 
oy = - Glt,x,y) (16) 


where 


Pitieoc7y) = ay 


and 
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G(t,xX,y) = bx 


giving 


>a 
ae ay Ci) 
d = a 
=r = — bx Ds 


with initial conditions 


C—O) eo =) sc (3) 


Ve 10) (4) 


II 
KK 


Also shown was the general form of the stochastic model: 


ee (t,m,n) 


ae Cee eee) + G(t,m,nt+l)P(t,m,nt1) 


— tenet) + 9G(t,m,n)) P(t, KR, nn) (20) 
with boundary conditions 
dP (t,m, Th ) 


—— 
aE PCE Th cel ha aE ei Ue e SDL ray 


+ ee pp bp (38) 


So 





= Ewe dB) Ss jase ara ) 


bp 


+ ei ane a +1) 


bp 
+ F(t m,n.) P(t,m,n,_) (39) 


np 


= ie ea Le) 
+ SUE OSL ag eas She UIE Se, 
= ea a (40) 
> 
n "yp 


ana initial conditions 


P(0,m,n) = 2 © (41) 
0 otherwise 


Substituting the square-law attrition rates gives the 


square-law stochastic attrition process: 


dP (t,m,n) 


aye —maiietcemit 7 thet ome ((t,m,nt+l) — (antbm)P(t,m,n) 


(19) 


and boundary conditions 
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op pg 
att SOE eo bp. + ee oo oes 
(41) 
adP(t,m,n,_) 
bp a 
a = an, P(t ,mrl yn) + SIE Ey Tig ahs Sree 
~ 2th off Wapalcih. | (42) 


| 


SE Ce ul ede) + Gh (Ree ae ey 


- Bol SOUT oP (43) 


ee initial conditions 


IL m = Mo and n= Dy 
Eto, n) = ¢ (44) 
|0 otherwise 


\ 


It was noted that the square-law comparisons have been 
made previously. However, as far as the author can tell, 
no previous comparisons have been made contrasting the models 
for breakpoint rules other than annihilation. Taylor [Ref. 30] 
has pointed out (and Adkins [Ref. 1] examined in further 
detail) that a neglected area in the modeling of land combat 
is the modeling of battle termination. In actual battles, 
several different events might cause battle termination. 
Bepone of the opposing forces is annihilated, then the battle 


ends; however, this is a rare event essentially never observed 
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in real combat, with a few exceptions (e.g., Iwo Jima, Alamo, 
etc.). If one of the opposing forces surrenders unconditionally, 
the battle will also terminate. The third and most common 

event which will result in battle termination is that of one 

of the opposing force breaking contact with the enemy and 
withdrawing from the battlefield [Ref. 13]. There are other 
possibie events which might result in battle termination, 

but they will not be considered. 

For the purposes of further discussion, a breakpoint is 
defined to be that state of a battle at which a unit con- 
Siders itself no longer capable of performing its mission 
and as a result elects to break contact with the enemy and 
withdraw from the battlefield. Therefore, when a unit with- 
draws from the battlefield strictly to avoid further combat, 
the unit is considered to have reached its breakpoint and 
lost the engagement. 

In many current land combat models a unit's breakpoint 
is determined by the percentage of casualties sustained by 
that unit [Ref. 1]. Although available empirical evidence 
shows that a deterministic battle termination model does not 
appear reasonable, it is an adequate "first cut" [Ref. 13]. 
And a breakpoint determined by the percent of casualties 
Sustained is more reasonable than one that assumes the battle 
goes to annihilation. For this reason, the investigation 
was done with deterministic breakpoints determined by the 


percent of casualties. 
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A. PROBABILITY OF WINNING 
Hypothesis 1-1 
As the breakpoint force levels, my Suge! Lgl are 
moved closer to the initial force Pievele. fi and 
n_, respectively; the difference in the 
probability of winning between the two models 
increases. 

Lee and Wannasilpa [Ref. 20] deduced that as the initial 
force levels are increased, the difference between the proba- 
bility of winning in the deterministic model and the stochas- 
tic model decreases. In fact it appeared that the difference 
disappears in the limit. However, the conclusion reached 
by Lee and Wannasilpa was based only on battles to annihila- 
tion. A look at the recursive relationship used in Appendix 
A to solve for the probability of winning is based on a 
summation from Hog 4 to mo Furthermore, the recursive 
relationships are from te (ae. ms These seem to indicate 
that the differences in the probability of winning decrease 
not just as m_ increases, but as m -™m,, increases. The 

O O p 
following figure seems to support this. Note, that for a 


deterministic model, the probability that Y wins, P has 


yo 
a Segoe nuoUs jump from one to zero as the ratio 
2 
(- cle P..) ) 
rah Bee through one. (Note: Fy and Fy are 


(1-(1-F,) 
the fractions of initial combatants each 


side is willing to lose.) For the square-law stochastic 
attrition process, however, Py moves continuously from one 


to zero reaching about 0.5 at a ratio of one. Note that as 


Py and FY decrease, Mp and Nop increase, so i hs and 
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ale decrease. As they decrease, Py force —StOChaStic 


model moves further and further from Py for the deterministic 


model. This was found to be true for all combinations of 


parameters considered. 


Hypothesis 1-2 


For fixed breakpoint casualty/initial force 
ratios (F,, and F,,), the difference between 

the probability Of winning for the deter- 
ministic model and the stochastic model 
decreases as the initial force levels increase. 


Lee and WannaSilpa [Ref. 20] deduced the above except 
their conclusion was restricted to battles to annihilation. 
The above hypothesis generalizes the conclusion to include 


any set of fixed Fy and Fy. Once hypothesis 1l-1l along with 


Lee and Wannasilpa's conclusion are accepted, hypothesis 
1-2 follows. The following figures support this. It is 


interesting to note that the plots of Py. are Similar for 


all fixed sets of FY and Fy. The only difference is in the 


angles of the curves. As Fy and F,, increase, the result 


ye 
is similar to rotating the curves counter-clockwise. This 
causes an increase in the difference between Py for the two 


models. 


Hypothesis 1-3 
As the forces move away from parity, the 
difference between the probability of 
winning for the two models becomes negligible. 


fits waesmenewitinel conclusion reached by Lee and Wanna- 


Silpa. Although their conclusion was for battles to 
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annihilation, the present findings do not alter the conclu- 
sion for other fixed force level breakpoints. The only 
change is in how far from parity the forces must be to make 
the difference negligible. 

The general consequence of this hypothesis is that unless 
the forces are near parity, there is little difference in 
the probability of winning between the two models. From 
this viewpoint, then, if one side is going to win decisively, 
there is little reason to use the more expensive stochastic 
model. Indeed, this has some intuitive appeal: The more 
decisively a side wins the battle, the less influence any 
uncertainties will have on the outcome. On the other hand, 
if the forces are near parity, the uncertainties in the 
evolution of the dynamics of combat may cause the expected 
loser to win. 

If the previous hypotheses are combined with this hypothe- 
sis, some interesting trends can be observed. It has been 


concluded that, with increasing Hh hs Shovel igual the 


o bp 
differences in the probability of winning between the two 
models decreases. When combined with hypothesis 1-3, this 
says that, with casualties on each side on the order of 20 
Or more, and with forces not near parity, there is no 
Significant difference between the probability of winning 
for each model. Even if the forces are near parity, if the 


force levels are very large, and both sides are willing to 


take 20 or more casualties; the difference in model win 
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probabilities is negligible. In these cases, if only the 
win probabilities were considered as measurement criteria 
between the two models, one would use the less expensive 
deterministic model as an acceptable Sei for the 
more realistic stochastic model. 

One must be cautious of this conclusion. It would be 
inaccurate to make this conclusion unless the only objective 
of the model is to determine which side will win. In most 
cases, there are other insights analysts try to gain from 
the model, such as how many casualties does each side take, 
at what points in the paceie are the most casualties being 
taken, how long will the battle last, etc. 

Nevertheless, the conclusions of this section do give 
indications of when the two models will not give the same 
Winners. In these cases, no further evaluation of the models 
is needed, because it is obvious the deterministic model is 
not an acceptable substitute. As Springall [Ref. 28] said; 

Probably the single most important 
comparison is between the two predic- 
tions of which side is going to be 
victorious. If the two methods cannot 
agree on this, there is little hope 
for agreement on the subsidiary 
attributes. 
feel ME HISTORY OF EXPECTED FORCE LEVELS 
If a determination has been made that the models will 


predict the same winner, further analysis may be made to 


see how well the models compare. Clark [Ref. 8] suggests 
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that a useful comparison is the time history of the expected 
number of survivors. In many models, one of the most utilized 
variables (performance and/or proxy variables) in measures 
of effectiveness for systems being evaluated by the models 
is the number of Survivors. It has been stated previously 
that a bias exists in the square-law and the linear-law. 
Lanchester himself suspected as much; he recognized that 
his differential equations were approximations to the casualty 
rates likely to be experienced in an actual battle. Quoting 
from Lanchester [Ref. 19]: 

Since the forces actually consist of a 

finite number of finite units (instead of 

an infinite number of infinitesimal units) 

the end of the curve must show discontinuity, 

and break off abruptly when the last man is 

reached; the law based on averages evidently 

does not hold rigidly when the numbers become 

small. 
Lanchester stated that his differential equations are based 
on averages implying an underlying stochastic process. He 
also suggests that his differential equations may be good 
approximations only so long as the force sizes are large. 


Although his reasoning is heuristic, Snow [Ref. 27] verified 


that the square-law is biased for small numbers. He showed 


that 
= = -an +a ££ nP(t,0,n) (37) 
n=1 
_ m 
dn —- 7 
ae = -bm + b Dee ne (em, Oo) (45) 
m=L1 
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The conclusion is that the solutions to Lanchester's equa- 
tions for the square-law are good approximations to expected 


Survivors so long as 


n 
e) 
ame nett,0,n) = 0 (46) 
n=1 
m 
fe) 
bee ome. (t ,m, 0) x0 (47) 
m=1 
Intuitively this seems sound. (46) and (47) express the 


boundary effects at time t. At early times in the battle, 

the expected force levels are near the initial force levels, 

and there is very little of the probability mass at the 

boundary. Thus the boundary effects are small. As the 

battle nears termination, though, much of the probability 

mass iS near the boundary and (46) and (47) will be signi- 

ficant enough to effect the process. As there are no boun- 

dary effects in the deterministic process, there will be 

a difference in force level trajectories of the two processes. 
Clark provided a derivation of (37) and (45). Additionally, 

he provided the following argument showing that the linear- 

law, which at first appears to be unbiased, is not. Through 


a Simple derivation, it can be shown for the linear-law: 


dm _ 
aaa aE (MN: t) (48) 
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= - bE(MN:t) (49) 


There are no bias terms in this expression (as opposed to 
(37) and (45)) since the random variable MN has a value of 
zero whenever M = 0 or N = 0. Consequently, the assumed 
differential equations for the linear-law appear unbiased. 
However, even though the analysis of the solutions to 
these equations show that the state solution is unbiased, 
the time solution interjects a source of bias. This is 


easily seen as a result of elementary probability theory: 


E(MN) = E(M): E(N) if an only if M and N 
are independent. 


In the linear-law, they are obviously dependent on each 
other. The rate of change of each force level 1s dependent 
on the product of the two force levels, causing the force 
levels to be dependent. 

Clark [Ref. 8] also gave numerical examples to show the 
amount of bias resulting from the difference in models. His 
results showed that there are significant amounts of bias 
resulting in both the linear-law and square-law cases. 

There are two things that should be noted about his results: 
First, his results are for small numbers of combatants; and 
second, his results are only for battles to annihilation. 


It was not shown in his dissertation that there is a 
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Significant bias with large numbers of combatants or with 
breakpoints other than annihilation. 

The following hypotheses are based on extension of his 
conclusions to include larger forces and fixed force level 
breakpoints. Appendix D is a derivation of the square-law 
stochastic attrition process showing that the bias noted 
by Snow exists with any set of fixed breakpoint force levels. 
Unfortunately, the form of the bias is not as clear as that 
for annihilation, since there are two additional terms. The 
use of the graphic aids noted earlier gives some indication 
that non-zero breakpoint force levels do not change the 
general progress of the battle. A look at the following 
figures will show that the time state probabilities progress 
in the same manner as in the case of battles to annihilation. 
Figures 12-15 are plots of the time state probabilities of 
the same model shown earlier except the breakpoint force 
levels are eight instead of zero (annihilation). Figures 
16-19 are the same except the breakpoint force levels are 
Sixteen. The differences are that the "wave" is absorbed 
by the boundary at Bina 2 and aeaat! instead of at (0,n) 
and (m,0). The question now is, What is the effect on 
the bias? 

A first analysis of (101) and (102) might give the 
indication that the bias is always positive (i.e., 
m(t)-x(t)> 0 for all t > 0). However, this is not true, 


as can be seen in the following argument. 
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Define 


eee) = c(t) (50) 
ey Cc) (51) 
then 
A, (t=0) = m(t=0)-x(t=0) = - any + an, = 0 C52) 
A, (t=0) = n(t=0)-y(t=0) = - bm, + bm. = 10 (53)) 
Let 
n 
S, (t) = at 2° EE pl “- By a ol Salt rd a a oy 
bp 
m 
re) 
+n PUG, 2 54 
bp 47? we (54) 
n 
m 
Oo 
S(t) = . LOE 6 (Mlg at Shale COM Sal mt 
np 
n 
e) 
+ m 2 PCE, yn) j (55) 
PP n=n +1 “bp 
bp 
af 
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<= = - ay (1) 
and 
am. an + S(t) (101) 
then 
= - cx = - an + SUA ay 
= Amex) =~ any) + sy (t) 
Substituting (50) 
dA, 
ae ot Sy (t) (56) 
and equivalently 
ae + S,.(t) (57) 
Gite X x 
with initial conditions 
(52) 
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Ay (t = 0) = 0 
Define 
Ay \ Smet) 
A= *| and [ee | z | 
\Ay j (Sy Ct) 
then 
A 0 a 
= - * U(t) 
b o|~ ~ 


and the solution of (58) is given by 


te 
A(t) =e ATA (t=0) we Rive eine dt 
s = 0 a 
where 
0 a 
i 
jon ?. 
Proof: 
A es tL = 
c= mien oeaon ee ur) = A fr e nte Tt) 
ii a 0 


il 


“AA(t) + U(t) 
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(53) 


(58) 


(59) 


Ua ar > AA (t) 





from Sylvester's Te Oren 


ane cosh (> ab t) -\fa/b sinn(\ ab t) 
e = 
-\;b/a sinh(>/ ab t) cosh (\/ab t) 
ieriCce 
t BR 
A, (t) = "4 S,(t)cosh(\/ab (t-t)) = Sy (t) \ a/b subeleie even (E15) |) ks 
(60) 
1 
at) = Pe S,(t)cosh( fab (t-t)) - S)(t)\/b/a sinh(jab (t-t)) at 
(61) 


If one now takes the case where one force is annihilated 


Pp = nD = 0), then Sy (t) and Sy (t) Simplify to 
n 
O 
S..(t) = © nP(t,0,n) (62) 
ie es 
n=0 
m 
O 
Sy (t) = Y mP(t,m,0) (63) 


which are simply the average number of Y survivors when X 
1s annihilated at time t and the average number of X sur- 


vVivors when Y is annihilated at time t respectively. It 


tcylvester's Theorem: For a polynomial function or 
convergent power series of a matrix, when the matrix has 
distinct eigenvalues, )., 


n nJ A-).1 
F(A) = 2 F(A;) 2 (=) where A iS an nxn matrix. 
L=] et eee ee 
j#1 
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seems intuitively appealing that if X wins decisively, 
Sy (t) + Q0. If this is true, then (60) will become negative 
and (61) will become positive. This seems to be confirmed 
by numerical results obtained. 

Another interesting case to examine is the case of 


complete parity (a=b, m_=n 


o Mo! We oy where Una ae” 


mm Chis case, Sy (t) = S, (t) = S(t) and 
f ; ome 
A(t) = «J SR keOsm ab (ts7))) — Sronc vicina (t- 1) ) dr 
(64) 
1 aioe shee 
Ay (t) = - See Teese ae (t>t)) > sinh(>/ab (t-1)) dz 
(65) 


By inspection, Ay (t) > 0 and Ay, (t) 2 Oe Or ate), and both 
biases are positive. This is reinforced by numerical 
results. 
The general conclusion reached from this analysis is 
that the biases are both positive unless one side is going 
to win very decisively. In this case, the bias on the 
Winner's side will be negative, while the loser's is positive. 
The remainder of this section will be a statement and 


discussion of hypcetheses concerning bias. 


Hypothesis 2-1 


Given fixed initial force levels and fixed 
attrition coefficients, as the breakpoint 
force levels increase, the numerical bias 
decreases. But aS a percentage of casualties 
of the deterministic model, the bias increases. 
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This hypothesis has heuristic appeal when the bias 


equations are analyzed. The bias terms are (see Appendix D): 


n 
O 
Sy (t) = ital BE Ep ey, + eRe ee 
bp 
™m 
O 
+n oy Pttem,n, _) } (54) 
DP m=m_ +1 ge 
np 
m 
O 


Sy (t) =i Jo) || 2 ee) + Lyrae Sia ay ous iia 


sot au 1 (55) 


Sete 


As ny and Mm, increase, the number of terms in each summation 
Pp Pp 
decreases. Although the second and third terms are multi- 


plied by n (which 1S increasing), it would seem that this 


bp 
may be more than offset by the decrease in number of terms. 
Even if this is true, aS a percentage, the bias will 
increase unless the actual bias decreases as fast as fs 
Although parametric analysis cannot logically prove the truth 
of this statement, it shows that the statement is certainly 
plausible. 
As all parameters except Mp and Nop are fixed, a 
numerical parametric analysis can be performed. A range 
of fixed parameters was chosen; small and medium values of 
Mo and Nos and several combinations of a and b. In all 


cases the results generally agreed. The following figures 


illustrate the results. They plot A,% and Aye as defined in 


x 


70 





100_% of Survivors on m Side of Total @asualties for Deterministic Battle 


90 

a=0.008 

b=0.004 

Mp=40 
80 Mpp=0 

No=40 

Npp=0 
70 
60 
50 
40 
30 

\ 
N 
\ 
20 ‘ \, 
x 
ow 
i \ : 
a 
‘ 


10 | \ 


: 2 w 6 460 410 | 


Percent of Time to Completion of Determini;' tc Battle 


Time Wistory Showing Rias 
FIGURE 20 


71. 


a 





100 4 of Survivors on m Side of Total Casualties for Deterministic Battle 


90 
=0.008 
M.=40 
{ Mop 8 
80 Nj=40 
70— 
| 
| 
| 
60 
50 
| 
40 
30 
, \ 
20 
10 
0-1 e pay See : = 
b 2b ad 6b ad 10b 
Percent of Tine to Completion of Deterministic 


Battle 


Time History Showing Pias 


PIGURE 21 


72 





100 


80 


70 


60 


50 


30 


20 


10 


% of Survivors on m Side of Total Casualties for Deterministic Battle 





20 Q 60 80 
Percent of Time to Completion of Deterministic . 
Battle 


Time History Showing Rias 


FIGURE 22 


73 


@ 





% of Survivors on m Side of Total Casualties for Deterministic Battle 


100 
90 
a=0 e 008 
b=0.0%%4 
60 AOS 
bp- 
N,=40 
Nbp=24 
70 
60 
50 
4O 
\ 
30 
20 
10 
0 
0 0 0 0 0 100 
Percent of Time to Completion of Deterministic 


Battle 


Time History Showing Pias 


FIGURE 23 


74 





100 % of Survivors on m Side of Total Casualties for Deterministic Battle 





D 20 4&0 60 80 100 
Percent of Time to Completion of Deterministic 
Battle 


Time Historv Showing Nias 


FIGURE 24 


7 





nis of Survivors on nm Side of Total Casualties for Deterministic Battle 


80 


70 


60 


50 


30 


20 


10 


= 


0 20 40 60 60 100 
Percent of Time to Completion of Deterministic Battle 





Time Vistory Showing Rias 


FIGURE 25 


iS 





10 


70 


60 


50 


30 


20 


10 


, % of Survivors on n Side of Total Casualties for Deterministic Battle 


SS 
SSM 
IN 


‘led ac seen SS ae - 
26 ub 60 86 106 : 


Percent of Tima to Completion of Deterministic 
Battle 


Time History Showing Bias 


FIGURE 26 


77 





* 


of Survivors on n Side of Total Casualties for Deterministic Battle 


100-4 
90 
=0.008 
b=0.00% 
M,=40 
80 Mp o= 16 
No=40 
Npp=16 
70 
60 
50 
40 
30 
20 \ 
10 
0 
; 20 4Q 60 fe 100 
Percent of Time to Completion of Deterministic 


Battle 


Time Wistory Showing Rias 


FIGURE 27 


78 





of Survivors on n Side of Total Casualties for Deterministic Battle 











100 


a 


a=0.008 
b=0.004 
M=40 
Mb p=24 
N,=40 
Npp=24 


80 
70 


60 


4O 


30 


0 8 100 


Percent of Time to Completion of Deterministic 
Battle 


Time History Sheowing Rias 


FIGURE 23 


qe 





10, 4 of survivors on n Side of Total Casualties for Deterministic Battle 


me N =h0O 


70 


60 


30 


20 





10 \ 


Percent of Time to Completion of Teterministic 
Battle 


Time Wistory Showing Pias 


FIGURE 29 


80 





a - eo omer = tm 


we eee ee eee 


A Survivors 
a2=02003 
b=0.004 
Mo=40 


wry: 
0 aac 


ae and a for Increasing Rreakpoints 


FIGURE 30 


ra ~ Survivors 
12 


24 3 


4: eft! Re Om Fei Vs iia son koa i 
A ' 





FIGURE 31 


81. 





Appendix C as functions of scaled time as defined in the 
Same appendix. The important things to notice are when the 
lines separate and how large the final separation is for 
the various input parameters listed on each figure. The 
figures are summarized in Figures 30 and 3l. 

In our examples, X is taken always to be the loser and 
Y to be the winner. As the breakpoints are increased, the 
difference in X Survivors decreases from slightly over three 
to slightly over one. This seems to confirm the first part 
of the hypothesis. Additionally, the rate of change of the 
bias is less than the rate of change of the breakpoints, so 
the percentage of bias increases, aS indicated in Figure 3l. 
This seems to confirm the second part of the hypothesis. 
However, the bias for the winning side does not follow the 
hypothesis (Y wins very decisively in this case). Indeed, 
it starts negative, becomes positive, and then decreases. 
But two things must be noted before the hypothesis is re- 
jected: First, the slope of the Ay curve 1s small; and 
second, it eventually decreases. This leads to the Ay3 
curve in Figure 31, where Ay® is increasing. 

It 1S argued in Appendix C that A% is the more important 
measure of the bias. If this 1s accepted, then the second 
half of the hypothesis is more important and, as indicated 
above, it is true. All of the results in Table 2 reinforce 
this. 

One final way of stating this conclusion is; the greater 


the differences m Tm cu Gierae 1) the less significant the 
O p O 


po. 


biases. 
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0.004 
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0.004 
0.004 
0.004 
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0.004 


0.004 
0.004 
0.004 
0.004 


0.004 


OS0015 
0.0015 


0.0015 
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0.0015 
0.0015 
Or0015 


050015 


M Mo 
40 40 
40 40 
40 40 
40 40 
40 40 
24 24 
24 24 
24 24 
24 24 
24 24 
40 40 
40 40 
40 40 
24 24 
24 24 
24 24 
24 24 
24 24 
TABLE II. 
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Hypothesis 2-2 
BOrmmEyxed FP. ana F., and fixed attrition 
rate coefficients, the larger the initial 
forces, the greater the numerical bias, 
but the smaller the percentage bias. 

A generally accepted conclusion for Lanchester-type 
models is a type of "law of large numbers". For a battle 
to annihilation, the differences between a deterministic 
model and the equivalent stochastic model "disappear" as 
the initial force levels increase. Hypothesis 2-2 is a more 
general statement of this not for battles to annihilation, 
but to fixed casualty/initial force level breakpoints. 
Additionally, it 1s concerned only with the differences in 
force level histories. Hypothesis 1-2 was a discussion of 
the "law of large numbers" for the probability of winning. 

There 1s an intuitive appeal for the hypothesis. The 
time history of the force levels for a deterministic model 
is something of an "average" time history of the stochastic 
model. As the number of combatants involved increases, the 
greater the aggregation, and the closer this average is to 
the mean force level history of the stochastic model. 

All but two parameters are held constant, so a numerical 
parametric analysis can easily be done. Figures 31 to 35 
show some of the results of the variation in parameters. 
Figures 37 and 38 summarize them and show that, for these 
fixed parameters, the hypothesis is true. 

Graphs of the results in Table III would give similar 


results and further reinforce the hypothesis. 
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The general result of this hypothesis is that, if the 
forces are of any appreciable size (greater than 20), with 
breakpoint force levels such that m —™m, and n_-n both 

O p Oo eDp 
exceed 20, the deterministic model is probably an adequate 


representation of the complex random process of combat. 


Hypothesis 2-3 
If any parameter is varied to bring the 
forces closer to parity, the time of the 
battle increases and the bias increases. 


This hypothesis is a result of an analysis of the solu- 


meron for Ay (t) and Ay (t) . As was shown previously, 


e Z poner 
A,(t) = “ [S.(t) cosh (>/ ab(t-1))-S, (1) \/ a/bsinh( yy ab(t-7) Jdt 
(60) 
c peer 
A(t) = a [S, (1) cosh ( \/ab(t-1)) -S, (1) \/a/bsinh (fab (t-1) ]dt 
(61) 


It appears that, as t increases, the integral will increase 
in absolute value. As only one parameter is varied, numeri- 
cal parametric analysis was easy. All results computed 
indicated the hypothesis is true. Figure 40 shows typical 
results. It 1s interesting to note that the hypothesis 
holds for the loser. However, it is not clear what happens 
for the winner, but it would appear that it is a continuous 


mamctiOon. if it is, then at some point, Ay reaches a 
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maximum negative value, and then increases to a positive 
value at parity. It then appears the hypothesis only holds 
for the loser. The results in Table IV show the same 


trends for the winner and loser. 


Hypothesis 2-4 
If a battle were to be terminated at a 
mimene-O.5t-, where t- 1s the completion 
time of the determiniStic model, there is 
no Significant difference between the 
results of the two models. 

Decisions to terminate battles are not always made 
because too many casualties have been taken. A good example 
of this would be the rear guard action of a division involved 
in a retrograde movement. The rear guard delays the enemy 
force, in many cases, long enough to allow the main force 
to retrograde and set up a defensive position, or at least 
allow a fresh rear guard to prepare for a further delaying 
action. In either case the rear guard commander will be 
tasked to delay the enemy for a fixed amount of time. When 
that time limit is reached, he will disengage, regardless 
of whether he feels he could engage the enemy longer. In 
fact, he may have taken very few casualties. This would be 
a Situation in which hypothesis 2-4 would be applicable. 

An evaluation of A, (t) and Ay (t) will show that they 
are generally increasing (in absolute terms), with increasing 
time. A study of various combinations of parameters has 


shown this to be the case. The question then is, at what 
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TABLE IV. 


tomes co bp x iy | fe 
24 19 24 0 0.88 0.35 27 
24 19 24 0 1.30 0.37 58 
24 9 24 0 1.68 0.29 91 
24 4 24 0 2.10 0.06 126 
24 0 24 0 2.48 -0.22 156 
6 0 6 0 1.06 -0.07 275 
We 0 6 0 3.22 1.93 493 
6 0 6 0 1.06 -0.07 275 
6 0 6 0 1.13 -0.10 291 
6 0 6 0 1.23 -0.14 312 
6 0 6 0 1.93 1.93 
6 0 2 0 1.06 -0.11 137 
6 0 1 0 1.18 -0.32 312 
6 0 12 0 1.34 -0.47 493 
6 0 12 0 ie S722 


Bias for Different Attrition Rate Coefficients 
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time into the battle do Ay (t) and A(t) become Significant? 
No method was found to determine this, but Figures 41 to 69 
are typical of all results found. 

As might be expected, the time at which the two time 
histories begin to diverge is strongly related to the final 
differences A,.3 (t-) and AVS (te) - As has been noted, as 
cor bp and Ws ays decrease, and the force levels are moved 
closer to parity, A,.3 (te) and Av% (t-) increase. The worst 
examples found did not diverge before 0.5 te, with one 
excaption. The one exception to this hypothesis is when 
the forces are at or very near parity. In this case, the 
hypothesis does not hold because, with the forces near or 
at parity, the rate at which casualties on both sides occur 
approaches zero as time increases. This gives, in a sense, 
an exaggerated time length of the battle and thereby distorts 


the relative amount of time A,% and A, % are appreciable. 


C. VARIANCE 

It is possible for the two models to agree on which side 
is going to win, and on the time histories of the average 
force levels, but the deterministic model may not adequately 
represent the complex random process of combat. This is 
because the outcome of the stochastic model is determined 
by the outcome of many uncertain events. The probability 
distribution function related to these uncertainties may 


cause the variance of the force levels to be considerable. 
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If the system being evaluated includes measures of effec- 
tiveness that are functions of the distribution of the final 
force levels, then in the case of large relative variance, 
the deterministic model is not an adequate representation 
of the complex random process of combat. 

Clark [Ref. 8] analyzed the variance of the square-law 
stochastic attrition process and hypothesized there are 
two types of time histories of the variance. The first is 
where the variance grows and converges asymptotically to a 
final level, while the second is where it grows to a peak 
level and then decreases asymptotically to a final level 
(see Figures 70 and 71). 

The qualitative pictures of the time state probabili- 
ties in Chapter II provide some useful insights into why 
the variance asymptotically approaches a value for every 
set of input parameters. As was noted, the change in the 
time state probability is a combination of diffusion and 
convective transport of the probability mass. The diffusion 
provides the variance in the force level. As the probability 
mass “hits" the boundary it remains there (is absorbed) ; 
for the probability mass absorbed, the diffusion effect 
has also ended and the variance of that probability mass 
is fixed. As more is absorbed, the less relative impact the 
remaining probability mass will have on the variance. Even- 
tually, as time grows, the remaining probability will have 


no effect, thus giving the asymptotic characteristic. 


i 2a) 


VARIANCE 
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However, because of the convective transport of probability 
mass, the remaining probability mass will not only continue 
to be diffused, it will, in general, arrive at the boundary 
closer to the breakpoint values (e.g., for the X force, the 
center of mass of the small probability mass, AP,, Baie a eaLinnes 
in the interval (t,ttAt,), will be at m survivors, but 


AP. Maniving Ln the interval (t+At t+At.) will be centered 


1! 
at m-k survivors). Thus, as time increases, the convective 
transport of the probability mass causes the average of the 
probability mass not absorbed to approach the breakpoint 
force level. The breakpoint force level acts as a boundary, 
and the probability mass cannot be diffused past that point. 
In fact, the probability that would have been diffused past 
the breakpoint "piles up" at it, and the variance is not 
guite as great as it might be. In fact, if the probability 
mass 1s transported fast enough toward the breakpoint, the 
variance of the remaining probability mass may decrease in 
Spite of the diffusion and might cause the overall variance 
to decrease. 

For illustration, take Figure 3 and collapse the 


probability mass to the X (m Survivor) axis. The result 


is: 
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This gives a qualitative picutre OL ee Ger uy so) m) ee the VarL— 

ance is calculated from the sum of the squared deviations 

from the mean. In this case, it appears the variance will be 

ierge” (in fact, it is on the order of magnitude of the mean). 
If we now take Figure 5 and collapse the probabl tasky 


mass to the X (m survivor) axis, the result.is: 


P(t=155.5,m) 








Pesci Ols |r? 10 0 


COUPAPSEO PRUBAILITY MASS 


PE yeiinenes 2 


Obviously the variance is much small than it was at 
t = 77.5. 

Note that X is the decisive loser in this battle; the 
convective transport has moved the probability mass toward 
the breakpoint fast enough to offset the diffusion effect, 
and the variance has decreased. On the other hand, if the 
same graphs were made for the Y force (the decisive winner) 
it would be obvious that the probability mass is not being 
transported toward the Y breakpoint very rapidly, and the 
variance continues to increase. 

Thus it seems that how the variance progresses is depen- 
dent on how rapidly the side is being attrited, relative to 
the other side. If a side wins, the variance will grow 
asymptotically, whereas if a side loses, the variance will 
peak and then decrease asymptotically. The results of all 
combinations of parameters in this study indicate this type 
of performance. 

Furthermore, the results in Table V indicated two 


general trends. 


Hypothesis 3-1 
is Mo Mbp increases, the variance increases 
eYic 


in num al value, but decreases in relative 
measure to the expected casualties. 


Hypothesis 3-2 


As the forces approach parity, the variance 
increases. 
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The second trend is in general agreement with previous 
conclusions in the respect that, unless the forces are near 
parity, the models are not significantly different. The 
first trend is in agreement with Hypothesis 2-1, so that 
the conclusions reached for it seem to apply also for the 
variance. Thus, even though the variance is a different 
measure for comparison of the two models, it seems to agree 


with previous conclusions. 
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IV. SUMMARY 


The intent of this analysis was to consider the question, 
can the complex random process of combat be adequately 
represented by a deterministic model. The models chosen 
for comparison were the deterministic and stochastic attri- 
tion versions of the Lanchester square-law attrition pro- 
cess. Three aspects of these processes have been compared 


and analyzed. 


Deterministic Model Stochastic Model 
A. Fixed winner and loser Probability of winning 
B. Time history of the Time history of the 
force levels expected force levels 
C. No variance of the Variance of the force 
force levels levels 


Conclusions in the form of hypotheses were discussed for 


each of the above comparisons: 


A. PROBABILITY OF WINNING 


Hypothesis 1-l 


As the breakpoint force levels, and 
an are moved closer to the initi 

force levels, m_ andn_, respectively, 
the difference fn the DEOoDaoie lL ty Or 
winning between the models increases. 


Hypothesis 1-2 


For fixed a ight casualty/initial force 
ratios (F,, and F,,), the difference between 
the probability ey winning for the deter- 
ministic model and the stochastic model 
decreases as the initial force levels increase. 
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Hypothesis 1-3 


As the forces move away from parity, the 
difference between the probability of 
winning for the two models becomes 
negligible. 


meen HESTORY OF EXPECTED FORCE LEVELS 


Hypothesis 2-1 


Given fixed initial force levels and fixed 
attrition coefficients, as the breakpoint 

force levels increase, the numerical bias 

decreases. But as a percentage of casual- 
ties of the deterministic model, the bias 

increases. 


Hypothesis 2-2 


For fixed F,, and F,,, and fixed attrition 
rate coefficients, the larger the initial 
forces, the greater the numerical bias, 
but the smaller the percentage bias. 


Hypothesis 2-3 
If any parameter is varied to bring the 
forces closer to parity, the time of the 
battle increases and the bias increases. 
Hypothesis 2-4 
If a battle were to be terminated at a time 
Meo, vnere &- 1S tie completion time 
of the eee nasitao model, there is no 


Significant difference between the results 
of the two models. 


Gee vARIANCE OF FORCE LEVELS 


Hypothesis 3-1 
As m_- increases, the variance increases 


in ntmeYvical value, but decreases in relative 
measure to the expected casualties. 
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Hypothesis 3-2 
As the forces approach parity, the variance 
increases. 
The hypotheses were formulated and evaluated through 
analyses of combinations of analytical solutions, numerical 


solutions, and graphical representations. 
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Ve. DISCUSSION 


The hypotheses formulated in this thesis provide a basis 
for the extension of several conclusions about the adequacy 
of the deterministic model. When one side is going to win 
decisively in the deterministic model, the models do not 
differ significantly and the deterministic model is to be 
preferred. How decisively one side must win to make the 
difference insignificant is a function of two sets of 
parameters; the initial force levels and the breakpoint 
force levels. As the initial force levels increase, the 
relative differences between models decrease as long as 
MMbp and No 2bp do not increase. Thus, with relatively 
large initial force levels (greater than 20) and breakpoint 
force levels such that there will be relatively large losses 
allowed (also greater than 20), the models are not signifi- 
cantly different. 

It must be recognized that these conclusions only apply 
to the very simple and idealistic Lanchester square-law 
attrition process and the square-law stochastic attrition 
process. As is always the case, the insights into combat 
dynamics obtained from these models are no more valid than 
the models themselves. However, intuition indicates the 
conclusions may be generalized. It is intuitively appealing 
to say that, no matter what the model or degree of complexity, 


if the deterministic model predicts that one side will win 
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decisively, the stochastic model will not disagree 
Sroniticantly. 

Nevertheless, it is felt that the hypotheses may 
provide a point of departure for the comparison of more 


complex models of combat. 
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VI. CONCLUSIONS 


HO@eecne 10dealized model {i.e., the so-called Lanchester 
“square-law" attrition process), it is concluded that the 
complex random process of combat can be adequately repre- 
sented by a deterministic model if the following conditions 
ese) MET : 

1. Each side starts with at least 20 combatants. 

2. Each side is willing to take at least 20 casualties. 

3. The forces are not near parity. 

4. If the forces are near parity, but each side 

initially has in excess of 40 combatants and 

is willing to take in excess of 20 casualties. 
Consequently, it seems plausible that similar results would 
be obtained if one were to compare analogous deterministic 
and stochastic “real world" combat models (i.e., large- 
scale complex system models of military forces). In other 
words, it seems reasonable to expect that if the opposing 
Military forces are not near parity and have relatively 
"large" numbers of initial combatants, a deterministic model 
may be used to model the dynamics of combat without losing 
any essential information about the battle's outcome. From 
the research reported here, the author concludes that one 
Nay safely take "large" numbers of initial forces (in par- 
ticular, of tanks), to mean units of battalion size or 
larger. Thus, it appears as though current U.S. Army deter- 


Ministic models such as the Bonder/IUA model are to be 
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preferred for sensitivity analysis because of their compu- 
tational advantages, once they are "calibrated" from Monte 
Carlo simulation output (for example, from DYNTACS or 


CARMONETTE) . 
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APPENDIX A 


COMPARISON OF TWO METHODS OF OBTAINING THE PROBABILITY OF WINNING 


In this appendix it is shown that two apparently differ- 
ent methods of obtaining the win probabilities are really 
the same. In his doctoral thesis, Springall [Ref. 28] 
derived a recursive expression for the probability of 
Winning in a homogeneous force, constant coefficient sto- 
chastic Lanchester model. His expressions and the develop- 
ment of them are quite different from those of Brown [Ref. 7]. 
Brown developed the familiar set of partial difference 


equations: 


A(r,s) 


Gee P(m,n:r-1,s) 


P(m,n:r,S) = 


Bir, s) 


aE A(r,s) + Blr,s) P(m,n:r,s-l) (66) 


with boundary conditions 


(67) 


oe 


0 otherwise 


1 m= r 
P(m,n, .:r/n,.,) = (68) 


0 otherwise 


Bests) = 0 n> Ss; 2 bot BBG, =QO mor 
(69) 
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fa =r andi n = Ss 
Pe sr s) = “bP Dp (70) 
“bp! bp 

0 otherwise 


where P(m,n mon) = probability X will win with m 


bp* 


survivors. When this is summed from thy to m,, it yields 


the probability X will win: 


Oo 
P. = 2 P(m,n 


2m 1) (ale) 
eon 


bps sc 


When the above equations are looked at carefully, what 
they are doing is explained by the following figure (Fig. 
74). 

The probability that X will win with m survivors is 
Simply the sum of the probabilities of all possible paths 


irom (m+n) eG) Vila Hs) Witte teem inOoadoblinktye On SOLNG  laem 


hye Es) 
A(r,s) + B(r,s) 


the probability of going from state (r,s) to state (r,s-l) 


BY .S) 
Piles) + Blr,s)- 


is simply the probability that X will win with Mm Survivors, 


state (r,s) to state (r-l,s) equal to and 


equal to And the probability X will win 
Or mo~1 SumvvOrs, OF <.., OF Mp survivors. 

A similar analysis of Springall's equations will show 
that the two methods are identical. Springall concludes 
that 
+1) (72) 


P(m,n 7M +N) = Bmp nt) xen 


bp bp 
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of Winning for Brown's Equations 


FIGURe 74 





where ICT a 


recursive relations; 


aia > aa. A(m_ ,n_) + B(m ,n_) 
o’ oO 0’ Oo 
B(mo,nt1)X(m_,n+1) 
O O 
A(m+1,n,)X(mt+1,n_) 
X{mno) = x(n.) + Btn) “np <™< % 
O O 
— A(m+1,n)X(mt1,n) + B(m,n+l) X(m,n+t1) 
X (m,n) az A(m,n) + B(m,n 
Define 
Fasc’, s) — \A(r,sS) + Bir;s))xX(r,s) 
Then 


I 
kr 


t 
P (mon) 


B(m_,n+1) 


' si ' 
Pp (mon) Atm, ntl) + B(m, ntl) Pp (m.,n+1) 


A(m+1,n_) 


' = ' 
E (m,n) A(m+1,n, see m+1,no = (mh1,n) 
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+1) is derived from the following set of 


(73) 


(74) 


i>) 


(76) 


a) 


(78) 


ie) 


(80) 





Paci nee (int | an) Bint Cin, nt) 
A(mt+tl1,n) + B(mt+l1,n) A(m,nt+l) + B(m,ntl) 


(81) 


P'(m,n) = 


Therefore 


AOE oe gu eS it +1) 


Pp Poe ; = 
(Mm oa no! A(m,,+1,n) + Bim, ti,n) (82) 
Aim, Nt 1)P (msn, tale) 
oe ote, Alin, +1) + 2B x Btn, FD LBS, 


What these equations say is, the probability of starting 


with (m..nQ) Survivors 1s one. The probability of having 


A(m,7n) 
one more casualty from X is and from 
Niner) = Bane nN) 
B(m en) Caso ©. 20 
Y is Ain, ze B(i yn) ° The probability of having (m,n) 


Survivors is the probability of having (m+l,n) survivors 
times the probability of an additional casualty being from 
X, added to the probability of having (m,nt+l) survivors 
times the probability of an additional casualty being from 
Y. This recursive relationship leads to the following 
figure (Figure 75), which is the same as Figure 74. Thus, 
even though the sets of relationships appear to differ, 
they do yield the same results. This was further supported 
by a comparison of results from Lee and Wannasilpa [Ref. 20] 
and the results from the enclosed comouter program, which 
was set up using the recursive equations of Springall. For 


the same parameters, the results were the same within 0.00001. 
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Method’ of Solving PeCuRSive Relationship for the Probability 
of Winning for Spvringall‘s Equations 


FIGURE 75 


The differences are assumed to have occurred because of 
roundoff errors within the computer system used. 


The program accuracy was further checked by 


™m n 
Oo O 


y: P(m ny 'M, n,Q) at. g 


P ( mM:m_,n_) = 1 (84) 
m=, eal n es “DP Oo. 06°08 
p 


Os 


All results were in the interval (0.9999692,1.0000000). 
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APPENDIX B 


NUMERICAL SOLUTION OF THE FORWARD KOLMOGOROV EQUATIONS 


Exact solutions for the time state probabilities for 
the square-law stochastic attrition process can be developed 
only for specific initial force levels. To the author's best 
knowledge, no solution for the time state probabilities for 
general initial force levels has been developed. Even if 
one were able to develop such a general solution, the result 
would be of no practical use due to its complexity. The 
solution for the linear-law stochastic attrition process 
was developed by Clark [Ref. 8], but it 1S too complex for 
practical use: So, solutions, if they could be developed, 
for more complex Lanchester-type stochastic attrition pro- 
cesses, would most likely be of no practical use. Therefore, 
the most practical way to solve such equations is through 
numerical integration techniques. The method chosen for 
this thesis was the fourth order Runga-Kutta method [Ref. 21]. 
The computer program is set up to solve Lanchester-type 


attrition process models of the form 


elie (t,m,n) 


SCCM eet BM, i))) bh Ce, ml, 2) 35) 
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with 


if) a No and m= m 
P(0,m,n) = (86) 
0 otherwise 


I 


One ec il (87) 


It 
ro) 
= 
A 
= 


P(t,m,n) 


and boundary conditions 


PS 
aP(t,m .,n) 
np 
= = A(m,t1,n)P(t,m, +1,n) (89) 


However, with a few minor alterations, the program will 
solve the more general Lanchester-type stochastic attrition 
processes with attrition rates of the form F(t,m,n) and 
eit, m,n) . 

Numerical precision of the Runga-Kutta methods is 
dependent on the length of the time step chosen. For the 
time steps chosen, numerical precision appeared to be 


satisfactory since 


Mm n 
O O 


» y P(t,m,n) was always in the interval 
m=Mm,, aise 


mom o992,1.0003). 
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Because the expected number of survivors was calculated 
from 


x O O 
if) ic) eo F yy mP (t,m,n) (90) 


m="bp nD bp 


it was felt the solutions were accurate enough. Additionally, 
because exact solutions exist for general initial force 

levels for P(t,m,n) and P(t,m,n.), they were used to 

check the accuracy of the approximate solutions. The 
approximate solutions never deviated from the exact solution 


by more than 0.00005. 
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APPENDIX C 


SCALING 


In a comparison of a stochastic model and the equiva- 
lent deterministic model, two classes of comparisons need 
to be considered. The first is the difference in results 
obtained from the two models. The second is to take these 
differences and see how they change with variations in 
parameters. One runs into difficulties when trying to make 
sense out of these comparisons when results are expressed 
in dimensional units (i.e., What is the significance of 
the fact that at 150 minutes into a battle, the difference 
in survivors is 1.8 tanks?). 

One method of overcoming this problem is the technique 
of scaling. If comparisons are made when results are in 
dimensionless quantities, it is much easier to see if there 
are any Significant differences. For example, one can 
eaSily compare two probability statements and note any 
difference between them. 

Scaling is a subset of a "science" called dimensional 
analysis, a technique often used in the physical sciences. 
For example, in the area of fluid dynamics, Reynolds deduced 
from some experiments that laminar flow broke down into 
turbulent flow at some critical velocity above that at 


Peel turbulent flow was restored to the laminar condition; 


itsee 





the former velocity being called the upper critical velocity, 
and the latter, the lower critical velocity. Reynolds was 
able to generalize his conclusions by the introduction of 

a dimensionless term called the Reynold's number. He found 
that certain critical values of the Reynold's number defined 
mrewupper and lower critical velocities for all fluids 
mrewing 1m all sizes of pipes, and thus deduced that single 
numbers define the limits of laminar and turbulent pipe 

flow for all fluids. Many other dimensionless numbers are 
used to generalize results — the Froude number, the Cauchy 
number, the Mach number, etc. Even in mathematics dimen- 
Sionless quantities are used to generalize results — 7 and 
e, for example. (7 is the ratio between the circumference 
of a circle and its diameter, regardless of the units of 
measure employed.) 

Scaling involves the ratio of two quantities with the 
Same dimensions. The dimension of a ratio of two quantities 
with the same dimensions is said to be dimensionless. The 
analysis in this paper 1s concerned with two classes of 
dimensions; time and quantity (number of casualties/survivors, 
and the differences in casualties/survivors). Either model 
(the stochastic or Pe termineistic) could have been chosen as 
the basis for forming ratios for the two classes of dimensions; 
the deterministic model was used as a basis. Time was trans- 
formed into a dimensionless quantity by using the ratio t/te, 


where te is the time at which the deterministic battle would 


Mey 





fooomeboeils CaSy to see that if 0 < t < te, then the ratio 
is between zero and one. 
The transformation of the quantities requires some 


explanation of the rationale behind it. The transformation 


used was 
A aes. number of casualties at t 
ratio number of casualties at ty for the deterministic mode 


It 1s easily seen that, for the deterministic model with 

< < . < < . 1 i 

Meo & < t OS ere ae Se A dimensionless quantity was 
desired that would express the "importance" of the number 
of casualties. A difference of one casualty out of forty 
is not as "important" as a difference of one casualty out 


erefrour. A expresses this "importance". For example, 
Pp p 


eat oO 


at te, eee = 0 for the deterministic model, while 

A patio x 0, in general, for the stochastic model. In fact, 
A - the difference in casualties at t, 
ratio 


number of casualties at ty for the deterministic model’ 


as the number of casualties at te for the deterministic 


model increases, A decreases, giving the deSired result. 


ratio 


Bes yS. 





ARE RPENDEX D 
DERIVATION OF EXPRESSION FOR BIAS IN THE SQUARE LAW STOCHASTIC 
Perret PON PROCESS WITH FIXED BREAKPOINT FORCE LEVELS 
From previous results it 1s known that the square-law 


stochastic attrition process is defined by: 


for oe <m < ™, and Nyptt en < nj) 
cote men) = anP(t,mtl,n) + bmP(t,m,ntl) - (antbm)P(t,m,n) 
(91) 
and on the boundaries of the state space, 
dP(t, mn.) 
aE = bmP (t,m,n +1) (92) 
ae (eym..,n) 
+ — = anP(t,m,+1,n) (93) 
also, 
P(t,m,n) = 0 fOr mo 7m and/or n>n) 
and 
1 m= Mm and n=n 
P(0O,m,n) = 4 
0 otherwise 
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Define 
Ae ftsmsn) = P(t,mtli,n) - P(t,m,n) 
and 
AP (t,m,n) = P(t,m,ntl) - P(t,m,n) 
so that (91) may be written as 


for mo otl <m MG and Dog Sie 


aP(t,m,n) 


ac = and P(t,m,n) = bmA_P(t,m,n) Vel) 


First consider two jointly distributed discrete random 
variables, M and N with probability mass denoted as 


Eu (men) - Further, assume that Eun Men) = Q for m < Mh 


and n < BD Then the expected value of M, denoted as 


E(M) is given by 


oO 


B(M) = 2 
eS if ap 


itm 8 
3 
Fh 
3 
Se 


This may be expressed in terms of the conditional expectation 


oS 


£m) 
aS E(M) = ¢f fy (n) E(M:N=n) 
n=n 
bp 
Note that 
m 
= O ne 
GMa. 5a 
an a rs = se mP(t,m,n) J (95) 
“bp bp 
dm : . 
EO compute a, which is given by equation (101) below, first 


sum equation (94) over alln, from Poy e to ny to obtain 


a No ares 
—— [ 2 P(t,m,n)] =a £ nA _P(t,m,n) 
dt n=n, +1 n=n, +1 a 
bp bp 
n 
O 
+ bm >» A_P(t,m,n) (96) 
n 
n=n, +1 
bp 
Note that 
n n 
O O 
r (eet m en) = [P(t,m,ntl) - P(t,m,n) ] 
m=n, +1 * n=n, +1 
bp bp 


P(t,m,n +1) = (ae oely ae, 


bp 


or 
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ve) 
7% Seeteem,n) = —- P(t,m,n, tl) 
1c lg ie eel a bP 
bp 
Thus 
n n 
d e) fe) 
——[({ f ECeeam,m) | = al & MAE tem.) | 
ans n=n, +1 N=noe tL ” 
bp bp 
= bmP (t,m,n, +1) (97) 
Adding (92) to (97) results in 
n 
re) n 
ad O 
——[{ 2 P(t,m,n)] = a[ 2 NACE Ce inn) | (98) 
age” m 
n=n. n=n, +1 
pp bp 


Multiply (98) by m and sum over m from ge to mo and 


add (93) multiplied by M5? 


WS ae - 
da _ dm 
: ak mP(t,m,n)] = ar 
“pp bp 
m n n 
Oo Oo Oo 
= + 
aa Pr a a a ee “+ ae = anP(t/T.,., ey) 
“DP bp bp 


(2) 


Observe that 


Uy ye 


a 1 


Oo 
iy mnd P(t,m,n) 
n=n, +1 
bp 
m n 
Oo Oo 
5 m, = ne (t ,mpw 
a i ah 
m a n 
Oo Oo Oo 
L 110 > nP(t,mtl,n) - £ ne(t;m,n)] 
aol We ve nS 
m n m n 
Oo Oo Oo Oo 
2 Meee nP(t,mtl,n) - £ m z ne Cem, a) 
I Pc 1 is n= ne! 
m +1 n 
Oo O 
» (m'-1) f gue Te Alin) pbek) 
m=M,' 2 SO rie 
as Mo 
= iy m e nr (t,m,n) 
— —- = —- 
m Mo Jk © api Mbp Jk 
We Seb n 
Oo Oo 
J (m*=-1L) SS Meee; ¥y) 
HE ole n=np ptt 
m n 
Oo Oo 
= vy InN 2 Met, ml, 0) 
m=m tt n=np ptt 
n 
fe) 
- zy Me (tem + Lm) 
“bP n=n, +1 
bp 
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Oo 
ae (m'-1) 2 MECe mm: ;n)} 
Baye n=n, +1 
m n 
O Oo 
- ue 2 me (t,m,n) 
Wt te n=n, +] 
n 
Oo 
os os >on aes” 
bp 
m n 
Oo Oo 
= Cie) TE ey) 
NE cs) n=ny ptt 
n 
Oo 
- P(t +1, 
Ms 28 Re ( MD n) 
bp 
Mo Ro 
= y nP (t,m,n) 
CNG th n=n ee 
n 
Oo 
- nP(t, +1,n 
Pp 
Substituting (100) into (99) 
— Mo are) 
= -a y > nP (t,m,n) 
oe m=m +1 n=n, +1 
“pp bp 
n 
Oo 
- ee z OE Sol eo 
n=ny ptt 
re 
+ am, ») met, me +.,n) 
Pp n=n, np 
Pp 


E59 


(100) 





™m nN 


O re) 
=-q f 5: nP(t,m,n) + amon Cee Sle Tle) 
ae +1 n=n, +1 ple “bP PP 
P bp 
m n m 
O O O 
==<a fF » UE Ge Tejon! de 5) OF ie ey eae, 
a _ bp bp 
= = = = + 
m Mh ln Np m Tp a 
+ SO sige Oya Ho sy 
m n n 
O O 
S oe y nP(t,m,n) + at nP(t,m, ,N) 
= n=n n=n P 
“pp ” “pp bp 
m 
O 
+at nh. PCE mn) 
2 bp bp 
m=m, t+ 
+ 
+ ees opt Ceo. 1/n,)) 
dm _ -an + alr nP (t nh) Eat tT 2) 
at a Mbp! bp bp Mbp Dp 
“bp 
™m 
O 
“i Np na ge ets at (101) 
p 
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Equivalently, 


m 
ai © 
ae bm + ce OES ont + eevee bp bp 
Pp 


2 2 PiGe, Payal (EO) 
“bP ke “bp 
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